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Abstract 

We explore the possibility that polaronic distortions in the paramagnetic 
phase of Lao.67Cao.33Mn03 manganites persist in the ferromagnetic phase by 
considering the observed electrical resistivity to arise from coexisting field- 
and temperature-dependent polaronic and band-electron fractions. We use 
an effective medium approach to compute the total resistivity of the two- 
component system, and find that a limit with low percolation treshold ex- 
plains the data rather well. To test the validity of this model, we apply it to 
the thermoelectric coefficient. We propose a plausible mean-field model that 
reproduces the essential features of a microscopic model and provides a com- 
parison with the experimental mixing fraction, as well as the magnetization 

and magnetic susceptibility. 
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There is a growing consensus that the so-called colossal magnetoresistivity (CMR) of 
La 1 _ x A x Mn03, where A is a divalent substituent, is larger when both double exchange 
and strong coupling to local lattice deformations are involved. |1|-§J In the double exchange 
model, [[U| electrons can hop from the singly occupied e2 9 orbital of Mn 3+ ions to empty 
&2g orbitals of neighboring Mn 4+ ions. Strong Hund's-rule coupling enhances the hopping 
matrix element when the S = 3/2 t2 g cores of the neighboring sites are aligned, thereby 
favoring ferromagnetism and an increased bandwidth. However, as Millis and coworkers 
miJ^] have emphasized, the Jahn- Teller effect in Mn 3+ , if strong enough, can lead to polaron 
formation and the possibility of self-trapping. The effective Jahn- Teller coupling constant 
X e ffi m this picture, must be determined self-consistently, both because it depends inversely 
on the bandwidth and because the effective transition temperature increases with decreas- 
ing X e ff. If X e ff is larger than a critical value A c , the system consists of polarons in the 
paramagnetic phase. As the temperature is lowered to the Curie temperature Tc, the onset 
of ferromagnetism increases the effective bandwidth, which reduces A e //, thereby increas- 
ing the effective transition temperature. As a result, the polarons may dissolve into band 
electrons if A e // drops below A c and the material reverts to a half- metallic, double exchange 
ferromagnet at low temperatures. 

In the Millis et al. model, the tendency toward polaron formation is monitored by a 
local displacement coordinate r, which is zero for X e ff < A c , and grows continuously as X e ff 
increases beyond that value. However, polarons are typically [§] bimodal-large or small-so 
that we should consider r to be a measure of the relative proportion of large polarons (band 
electrons for which r f=a 0) and small polarons (for which r is an atomic scale length). In- 
deed, there is growing experimental evidence f9HT2|| that polaronic distortions, evident in 
the paramagnetic state, persist over some temperature range in the ferromagnetic phase. 
This paper explores that possibility, in the effective medium approximation, by considering 
the observed electrical resistivity to arise from a field- and temperature-dependent concen- 
tration c(H, T) of metallic regions within a semiconductive, polaronic background having 
an activated electrical conductivity. |L3| We test the validity of this model by applying it 
to the thermoelectric coefficient, also within the effective medium approximation, and ex- 
tract c(H, T) for stripe-like domains. We propose a mean-field model, in which c(H, T) is 
a secondary order parameter, that reproduces the qualitative features of the experimental 
data. 

The La 2 /3Cai/ 3 Mn03 film samples used in this study were prepared by pulsed laser de- 
position onto LaA10 3 substrates to a thickness of 0.6 fim. As described previously |TJ|], they 
were annealed at 1000 °C for 48 hr. in flowing oxygen. Measurements were carried out in 
a 7 T Quantum Design Magnetic Property Measurement System with and without an oven 
option provided by the manufacturer. A modified sample rod brought electrical leads and 
type-E thermocouples to the sample stage. A bifilar coil of 12 fim Pt wire was calibrated 
to serve both as a thermometer and to provide a small heat input for the thermopower 
measurements. Measurements in fields up to 7 T could be carried out over the tempera- 
ture range 4 K < T < 500 K. Following the transport measurements, magnetization data 
M(H, T) were acquired up to 380 K by conventional methods. 

Figure 1 shows the magnetization and resistivity data over the temperature range of 
interest. The data below 200 K exhibit essentially field independent, metallic behavior, and 
are well fit by a power law, 
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pu(T) = [0.22 + 2 x 10"° KT 2 T 2 + 1.2 x 10~ 12 RT 5 T 5 ] mtt cm 



(1) 



Above 260 K, the resistivity is exponential and again independent of field, given [15| by 
the form expected for the adiabatic hopping of small polarons, 



p M {T) = (lApQ cm K _1 )Texp(1276 K/T). 



(2) 



Both limits are displayed in Figure 1 as dotted lines. Our assumption is that these represent 
the resistivity of band electrons and polarons, respectively, and that the transition region 



can be represented by the conductivity of an effective medium [IB] characterized by a mixing 
factor c(H, T) which we envisage to be the fraction of the carriers that are in the metallic 
state. We assume the effective resistivity p(H, T) to satisfy the general expression for ellip- 
soidal metallic regions (resistivity pn(T)) embedded in a semiconducting matrix (resistivity 
PhtiT)) H, 



c - Pit) 



.Pit + 9\\(p ~ Pit) Pit + 9±(p - Pit) 
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(3) 



Here, g± ~ 1/2, g\\ = (b 2 pht/o?p) ln[l + (ap/bpht)} and a, b are the (major, minor) axes of the 
prolate ellipsoids. As a first approximation, we set c(0,T) = M(0,T)/M sat , using the data 
of Fig. la. The chain curves in Fig. lb show that the magnetization roughly reproduces the 
shape of the resistivity curves at zero field. That the metallic resistivity switches on quickly 
as the system orders indicates a low percolation threshold and a large aspect ratio (a >> 6). 
This suggests that the metallic regions appear as stripe domains, rather than as compact 
clusters. The magnetization does not approximate, at all, the relative concentrations in 
applied fields, as seen in the solid line in Fig. lb. As an alternative, we solve Eq. (3) for 
c(H,T) using the experimental resistivity p exp (H,T) , with the result shown in Fig. 2a for 
a/b = 50. 

As a second test of this approach, we consider the Seebeck coefficient S(H, T), measured 
over the same temperature range, and plotted in Fig. 2b. We fit the low temperature 
thermopower arbitrarily to a power law 



S lt (T) = [(0.051 K -1 )T - (1.3 x 10~ 4 RT 2 )T 2 - (3.2 x 10~ 7 K" d )T d ] pV /K 



and the high temperature data |TjJ to the form expected for small polarons, 

S ht (T) = [(9730 K)T- X - 29] pV/K. 



(4) 



(5) 



Low and high temperature limits are displayed as dotted lines in Figure 2b. A remarkable 
feature of the transport properties of composites is that knowing any two effective transport 



properties and those of the constituents determines the third. [0 A closed expression for 
the Seebeck coefficient has been derived by Bergman and Stroud |18|] which, by recognizing 
that the effective conductivities match the high and low temperature values in the respective 
limits, we can write in the following form, 



S(H, T) 



(Su - S ht ) (Kg - ft/tt) Pit Pht 
(Kht Pht - Kit Pit) (Pht - Pit) 
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(6) 
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Here, Ku, km, and k(H, T) are the thermal conductivities of low and high temperature phases 
and the effective medium, respectively. While small changes in k(H, T), of order 20% in the 
transition region, have been reported [HJ , the thermal conductivity of this thin-film sample 



is dominated by the substrate. All three thermal conductivities are, therefore, equal at each 
temperature-field point, enabling us to simplify the expression to, 

S (H, T) = [phtSu - puSht + Pcxp (H, T) (S ht - Sit)]. (7) 

Pht ~ Pit 

The results are shown as solid lines in Fig. 2b; the agreement is excellent, with no adjustable 
parameters. Measurements using single-crystal samples are underway to explore the effect 
of k(H, T) on the analysis. 

The precise shape of the c(H, T) curves depends on the form of the effective medium 
model chosen, being sharper the lower the percolation threshold. Nonetheless, it is clear 
that the metallic fraction increases rapidly below the transition temperature, but does not 
broaden in applied field to the same extent as the magnetization. It is useful, therefore, to 
explore under what circumstances this can occur. The essential feature of the Millis et al. 
model is that the effective Jahn- Teller coupling constant X e ff is very near its critical value 
in the paramagnetic phase. There are two phase transitions to consider a zero-temperature 
polaronic transition at a critical coupling constant A c that would occur in the absence of 
Hund's rule coupling, and a double-exchange ferromagnetic transition at Tc for A e // = 0. 
The dependence of the bandwidth on magnetic order in the double exchange picture causes 
both X e ff and T c to depend on the temperature through the magnetization. We propose 
here a simplified mean-field model that reproduces the essential features of the microscopic 
calculation to demonstrate that c(H, T) is a secondary order parameter, driven by the pri- 
mary order parameter, namely the magnetization, and that it resembles the experimentally 
determined curves. We assume that the ferromagnetic free-energy functional is of conven- 
tional form 

F rnag = ^(T/T c -l)m 2 + ~bm*-mh, (8) 



written in units of |2fJ 3Sk B T c /(S + 1) = lMk B T c for S = 2(1 - x) + 3x/2 = 1.83 
and x = 1/3. Further, we have h = gps(S + l)H/3k B T c = -ff/2360 kOe. The secondary 
parameter c(H,T) is asumed to be driven by the difference X e ff — A c . We approximate the 
dependence of X e ff on the magnetization by writing X e ff — A c a a — 7m 2 + where a is 
small and positive. The electronic free energy can then be written, in the same dimensionless 
units as Eq. (R), as 

F el = l -{a-W)S + \pc\ (9) 

Minimizing the total free energy F = F mag + F eh we obtain 

dF 

— = (T/Tc - 1 - 7c 2 )m + bm 3 - h = (10) 
dm 

OF 

— = ( a - ^m 2 )c + (3c 3 = (11) 
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From Eq. (|TTD it is obvious that the concentration of metallic electrons is zero until the 
magnetization reaches the value m = *J~otpy, beyond which point c increases. In the limit 
a — > 0, c is proportional to m, with the result that b — > b — 7 2 / (3, signalling a tendency for the 
system to approach a tricritical point and first order transitions as the coupling constant is 
increased. Note that the existence of a non-zero concentration c can be considered to increase 
the critical point to (1 + 7c 2 )Tc, causing the magnetization to increase more rapidly than 
would be the case without coupling to the metallic electron concentration. 

To proceed, we recognize that Eq.fllOD is the expansion of a Brillouin function, 

m = B s + 7C 2 )m + h)j , (12) 

while the equation for c can be assumed to be a small-c expansion of 

c = tanh [(1 — a + 7m 2 )c] . (13) 

In Fig. 3 we show the simultaneous solutions of Eqs. ([12] & 0) for a = 0.02 and 7 = 0.3 
at H = 0, 24 kOe, and 48 kOe. Application of the magnetic field increases the temperature 
at which c becomes non-zero by 7% or 20 K, consistent with the experimental data in Fig. 
3 and, unlike the magnetization, does produce a high-temperature tail in c. As no thermal 
factors are included in the definition of c, the concentration of free carriers does not approach 
unity, and therefore differs slightly from the experimental c(H,T) obtained from Eq. @. 
The abrupt appearance of band electrons in this model produces a kink in the zero field 
magnetization curve at the onset temperature, seen as a deviation from the H = 0, 7 = 
curve. The inset to Fig. 4 shows the inverse of magnetization data in an applied field of 1 kG 
taken on a single crystal of Lao.7Cao.3MnO"3, |21| along with H/m for the 24 kG data of Fig. 



3. The experimental kink moves to higher temperature as the field is increased, as predicted 
by the model, but occurs further from the actual transition than the model predicts. The 
transition temperature for this sample, obtained by a scaling analysis of the data below the 
kink temperature, is 216 K. Clearly, the mean-field model proposed here cannot capture the 
precursive behavior due to critical fluctuations, important in the actual experiment, that 
would move the kink to higher temperatures. 

As noted above, a number of other experimental probes have suggested the coexistence 



of localized and delocalized d-holes in the ferromagnetic state. Booth et al. ['I L] have defined 
the density of delocalized holes n^h analogously to our Eq. (|^) from the width of the Mn-0 
EXAFS peak. They suggest that n dh oc exp(3.5m), which leaves a finite density of delocal- 
ized carriers above Tq and a nearly linear increase in n^h below. The fits for the mean-square 
width of the Mn-0 bond length distribution, however, look remarkably like the magnetiza- 
tion itself. Louca, et al. 0], in a Sr-doped sample with comparable Tc, report a gradual 
change in the number of nearest Mn-0 neighbors as the temperature is reduced below the 
critical point. They interpret their data in terms of a gradual transition from single-site to 
triple-site polarons, rather than to free carriers as favored here. Similar neutron data by 
Billinge, et al. show that the 0-0 bond-length distribution begins to narrow below Tc 



in a manner similar to our c(0,T). Recent Raman results [T(| find two contributions in the 



Raman intensity a diffusive component associated with small polarons and a continuum con- 
tribution from free carriers. The relative intensities of these components have temperature 
dependences that are very similar to 1 — c(0, T) and c(0,T), respectively. 
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That the field and temperature dependence of the thermopower can be deduced directly 
from resistivities and field-independent component Seebeck coefficients lends strong support 
to a picture in which conducting regions arise and percolate in the vicinity of the phase 
transition. The rapid decrease in resistance requires a low percolation threshold which in 
turn, points to stick-like, or stripe-like domains. Our simple mean-field model ignores 
a number of features that should be included in a complete treatment. In particular, we 
have excluded a term m 2 c because it leads to a first-order transition for all values of the 
parameters; we cannot rule it out on symmetry grounds. Similarly, there should be a mix- 
ing entropy in the electronic free energy which, at sufficiently high temperatures, will lead 
to thermal dissociation of the polarons. Finally, we have not included gradient terms and 
therefore ignore inhomogeneous thermal fluctuations that are certain to be significant in a 
system such as this where there are competing order parameters. Nonetheless, our phe- 
nomenological approach provides a qualitative understanding of the field and temperature 
dependence of the transport properties while correctly predicting the existence of kinks in 
the magnetization curves. 

This work was supported in part by DOE Grant No. DEFG0291ER45439 through the 
Illinois Materials Research Laboratory. MJ acknowledges support from U.S. Department of 
Energy at Los Alamos National Laboratory, NM. 
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FIGURES 



FIG. 1. (a) Magnetization data for this sample, (b) Resisitivity data as functions of field 
and temperature. The broken curves are effective field calculations at two values of a/b, assuming 
that the low field magnetization represents the metallic concentration. The solid line uses the 
magnetization measured at 10 kOe, and a/b = 50. Dotted lines are the low and high temperature 
limits discussed in the text. 

FIG. 2. (a) The mixing factor c(H, T) extracted from the resisitivity using the effective medium 
approximation with a/b = 50. (b) Seebeck coefficient data and results of a computation using the 
measured resistivity p exp (H, T) as described in the text. Dotted lines show the low and high 
temperature fits. 

FIG. 3. (a) The mixing factor c(H, T) calculated in the mean-field model with a = 0.02 and 
7 = 0.3. (b) The magnetization calculated with the same parameters. The dotted line shows the 
non-interactive case for comparison. 

FIG. 4. Inverse magnetic susceptibility H/m vs temperature near the M-I transition for H = 24 
kOe. The dashed line is the zero field extrapolated behavior. The appearance of free carriers induce 
the rise of the effective Tp, leading to a kink in the susceptibility. Inset Inverse susceptibility of a 
single crystal sample of Lao.zCao.sMnOs, measured in a field of 1 kG. 
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